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Abstract 

Using the recently introduced boundary form factor bootstrap equations, we map 
the complete space of their solutions for the boundary version of the scaling Lee- Yang 
model and sinh-Gordon theory. We show that the complete space of solutions, graded 
by the ultraviolet behaviour of the form factors can be brought into correspondence 
with the spectrum of local boundary operators expected from boundary conformal 
field theory, which is a major evidence for the correctness of the boundary form factor 
bootstrap framework. 

1 Introduction 

The bootstrap program aims to classify and explicitly solve 1+1 dimensional integrable 
quantum field theories by constructing all of their Wightman functions. For bulk theories, 
the first stage is the S-matrix bootstrap: the scattering matrix, connecting asymptotic in 
and out states, is determined from its properties such as factorizability, unitarity, cross- 
ing symmetry and the Yang-Baxter equation supplemented by the maximal analyticity 
assumption which results in the complete on-shell solution of the theory, i.e. the spectrum 
of excitations and their scattering amplitudes (for reviews see [HE]). The second step is 
the form factor bootstrap, which allows one to determine matrix elements of local oper- 
ators between asymptotic states (form factors) using their analytic properties originating 
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from the already known S-matrix. The form factors are then used to build the correla- 
tion (Wightman) functions via their spectral representations, yielding a complete off-shell 
description of the theory (see [3] for a review). 

The first step of an analogous bootstrap program for 1+1 dimensional integrable bound- 
ary quantum field theories, the boundary R-matrix bootstrap, was developed in the pio- 
neering work by Ghoshal and Zamolodchikov [4j; it makes possible the determination of 
the reflection matrices and provides complete description of the theory on the mass shell. 

For the second step matrix elements of local operators between asymptotic states have 
to be computed. In a boundary quantum field theory there are two types of operators, the 
bulk and the boundary operators, where their names indicate their localization point. The 
boundary bootstrap program, namely the boundary form factor program for calculating 
the matrix elements of local boundary operators between asymptotic states was initiated 
in [50. 

In this work we make further progress in understanding the boundary form factor 
bootstrap by classifying and counting the solutions. In the bulk case it was proposed by 
Cardy and Mussardo that the space of solutions of form factor axioms is to be identified 
with the space of local operators |8j. The program of counting the solutions was performed 
for several interesting models p2 [10], HH H2] and it was shown that the results agree with 
the spectrum of local operators expected from the Lagrangian description in terms of the 
ultraviolet limiting conformal field theory. Here we take the models treated in [5] and 
perform a similar counting of solutions, which can then be brought into correspondence 
with the spectrum of boundary operators predicted by boundary conformal field theory. 

The outline of the paper is as follows. In Section 2 we give a short overview of the 
boundary form factor axioms and also introduce some basic notions concerning form factor 
solutions. In Section 3 we investigate the scaling Lee- Yang model with boundary, which 
has two physically inequivalent boundary conditions. We start with the simpler boundary 
condition (denoted by 1) and develop all the concepts necessary to perform the counting of 
the form factor solutions. Then we apply these to the other boundary condition (denoted 
by $) too, while in Section 4 we treat the boundary sinh-Gordon model using the same 
tools. We give our conclusions and discuss major open problems in Section 5. 

2 The boundary form factor bootstrap and the classifi- 
cation of solutions 

2.1 The boundary form factor axioms 

The axioms satisfied by the form factors of a local boundary operator were derived in [5]. 
Here we only list them without much further explanation. Let us suppose that we treat 
an integrable boundary quantum field theory in the domain x < 0, with a single scalar 

1 There exists no analogous framework for bulk operators in the presence of the boundary yet. Large 
distance expansion for their correlators in terms of bulk form factors can be given using the boundary 
state formalism, but so far this approach has only been used to compute one-point functions (SIEJ- 
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particle of mass m, which has a two-particle S matrix S(B) (using the standard rapidity 
parametrization) and a one-particle reflection factor R{9) off the boundary, satisfying the 
boundary reflection factor bootstrap conditions of Ghoshal and Zamolodchikov [4]. For a 
local operator 0{t) localized at the boundary (located at x — 0, and parameterized by the 
time coordinate t) the form factors are defined as 

out (#1 j 6*25 ••• 5 ^ml^Wl^l) ^2, • • • , 0n)in = 

F° (f)' B' B' ■ 0i Bn B )p- irnt (52 cosh 6i-22 cosh &'.) 

for &i > B 2 > ■ ■ ■ > B n > and B[ < B' 2 < . . . < 9 m < 0, using the asymptotic in/ out 
state formalism introduced in [13]. They can be extended analytically to other values of 
rapidities. With the help of the crossing relations derived in [5] all form factors can be 
expressed in terms of the elementary form factors 

out(O|O(O)|0i, B 2 , ... , B n ) in = F n (9i, 02> • • • > n ) 

which can be shown to satisfy the following axioms: 

I. Permutation: 

F n (0i, • • • , 0j, 0i+i, ■ ■ ■ , n ) = S(Bi - 9 i+ i)F n (9 1 , . . . , 6 i+1 , 9i, . . . , 9 n ) 

II. Reflection: 

F n (9x, . . . , n _i, B n ) = R(B n )F° {B\, . . . , 0„_i, — n ) 

III. Crossing reflection: 

F°(9 h 2 , ... , n ) = R(tK - 9 1 )F%{2m -9 h 9 2 ,..., 9 n ) 

IV. Kinematical singularity 




or equivalently described as 




V. Boundary kinematical singularity 



-i Res F° +1 (9 + ^,B 1 ,...,B n ) = 9 -{l-H S(j - 9^^ ...,9 n ) 
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where g is the one-particle coupling to the boundary 

VI. Bulk dynamical singularity 

-i ResFf +2 (0 + iu, B' - iu, 9,,..., 9 n ) = TF° +1 (9, 1; . . . , 9 n ) 
e=e 

corresponding to a bound state pole of the S matrix 

(in a theory with a single particle, the only possible value is u = 7r/3). 

VII. Boundary dynamical singularity 

-i Res F n o +1 (0 x , 9 n , 9) = ~gF°{9 x , 9 n ) 

9=iv 

which corresponds to a pole in the reflection factor describing a boundary excited state: 

R(0)„?tll , 9~iv 
9 — iv 

We further assume maximum analyticity i.e. that the form factors have only the minimal 
singularity structure consistent with the above axioms. 

2.2 Solution of the axioms 
2.2.1 The general Ansatz 

The general form factor solution can be written in the following form [5] 

n 

F n (9t, 9 2 ,..., 9 n ) = G n (0!, 0* ... , 9 n ) J] r(9i) J] /(0, - 0,)/(0 i + 9 3 ) (2.2) 

i=l i<j 

where / is the minimal bulk two-particle form factor (2PFF) satisfying the conditions 

f(9) = S(9)f(-9), f(iir + 0) = f(iK-e) 

and having the minimum possible number of singularities in the physical strip < 9 < tx 
together with the slowest possible growth at infinity [H], and r is the minimal boundary 
one-particle form factor (1PFF) satisfying 

r(9) = R{9)r{-9) ; r{m + 0) = R{-9)r{m - 9) (2.3) 

plus analytic conditions similar to those of /, but in this case in the strip < 9 < ir/2. 
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The functions G n are totally symmetric and meromorphic in the rapidities They 
are also even and periodic in them with the period 2ni, so they can only be functions of 
the variables 

y t = e e ' + e~ Bl 

Let us now turn to the analysis of the singularity structure. In a theory with only one 
particle, the only possible singularity of the S matrix in the physical strip is located at 
6 = 2ni/3 corresponding to the self- fusion of the particle (plus the crossed channel pole 
for the same process at ni/3) and the relevant fusion coupling is defined as 

T 2 = -i Res S{9) 

We suppose that the 2PFF function / is chosen such that it has a pole at 6 = 2ni/3 so 
that it encodes this singularity (for an example see (13.21) ). We further assume that the 
boundary dynamical singularities (but not the kinematical ones!) are similarly contained 
in the 1PFF function r. Then the functions G n only have singularities at the positions of 
the kinematical singularities, and so they can be written in the form 

n la a a \ Qn(VuV2-- .,y n ) 

(*n(Vl, t>2, ■ ■ ■ , V n ) ~~ 



UiVi Yl(vi + Vj) 



where the Q n are entire functions symmetric in their arguments. 

2.2.2 Two-point functions and scaling weights 

The two-point functions can be written using a spectral representation: 

(0|C(t)C(0)|0) = V — — / d9 x de 2 . . . d6 n F n F+ exp -tmtY] cosh^ 



where time translation invariance was used and the form factors F ni are defined as 



(2.4) 

a factors F n , are defined ; 

F n = (010(0)1^, e 2 , ... , e n ) m = F^(e x , e 2 ,..., e n ) 

and 

F+ = in (6i, 9 2 , ... , 0„|0(O)|O) = F®{iTT + 9 n , m + n _ 1; ...,m + 9 x ) 

which, for unitary theories, is the complex conjugate of the previous one: F+ = F*. In 
the Euclidean time r — it the form factor expansion of the correlator converges rapidly for 
large separations since multi-particle terms are exponentially suppressed. 

We are interested in operators which can be classified according to their scaling di- 
mensions, which means that the two-point function must have a power-like short-distance 
singularity 

(0|O(r)O(0)|0> = ^ (2.5) 
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where A is an exponent determined by the ultraviolet scaling weights of the local fields. 
More precisely let us consider the conformal operator product expansion 



hi 

where h is the ultraviolet weight of the field O and the hi are the weight of the Oi. It is 
obvious that 

^ \ — — // *^Tnin 

where h m { n is the minimum of the weights of the operators appearing in the expansion. In 
many cases it is h m[n = (corresponding to the weight of the identity) and then h and A 
are identical. 

It is well-known that the presence of a power-like singularity means that the form factors 
themselves can only grow exponentially in the rapidity variables [15] and so the functions 
Q n are symmetric polynomials in the variables yi. The singularity axioms give recursion 
relations for the polynomials Q m for which we introduce the abbreviated notation 

Qn — [Qn+2\ 
Qn = D [Qn+l] 

Q n = B[Q n+1 ] (2.6) 

where /C, T> and B denote the recursion relation resulting from the bulk kinematical, bulk 
dynamical and boundary kinematical singularity axioms respectively. We give explicit 
forms of these relations for specific models later. 
If the form factors grow asymptotically as 

\F n (9 l + e,9 2 + e,...,9 n + 9)\ ~ e xW as |6»| ^ oo (2.7) 

(where we assume that the exponent x is independent of the level nj^ then the leading 
short-distance behaviour of the individual terms can be easily shown to be 1/t 2x and so 
naively A = x. In fact, leading logarithmic corrections in r can (and in many cases do) 
sum up to an anomalous contribution to the ultraviolet exponent A and therefore we call 
x the naive (or engineering) dimension of the form factor solution. 

2.2.3 Towers and local operators 

Following the ideas in the work of Koubek and Mussardo on classification of the bulk form 
factor solutions [9], the set of elementary form factors corresponding to a local operator 
forms a tower of form factors, graded by the number of particles: 

2 For all models considered here we show later that the space of solutions is spanned by so-called 'simple 
towers' for which the asymptotic exponent x is indeed independent of the level n. 
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Due to the recursion relations (12.61) the solutions (12.21) 



Qn(yi,y2 ■■■,y n ) 
YliVi Yl(yi + yj) 



n 



F n (9l, 02, ■ ■ ■ , n ) 



i<j 



i=l i<j 



of the form factor equations can also be classified into towers which consist of a single form 
factor at each level n such that the solutions at different levels are linked together via the 
relations (12.61) . As we show later in the explicit examples, the recursion relations are such 
that the naive scaling dimension x is independent of the level n and therefore it can be 
assigned to the tower itself. Therefore the space of form factor solutions can be rearranged 
into a space of towers T, which can in turn be graded by the naive scaling dimension. Due 
to the linearity of the form factor axioms, the space of towers is a linear space i.e. whenever 
T = {F n } ngN and T' = {F n } nGN are two towers, so is their general linear combination 



While the space of towers (T) is always infinite dimensional, the subspaces corresponding 
to a given value of x can be finite dimensional; in that case it makes sense to count the 
towers graded by x. We can introduce the character of this space by 



where d x is the linear dimension of those towers which have naive dimension x. 

The space £ of scaling local boundary operators is also a graded linear space, where 
the grading is given by the ultraviolet scaling weight A. The character of that space 
can be obtained from boundary conformal field theory. Our aim here is to bring the two 
graded spaces into correspondence with each other via comparing their characters, which 
is actually nothing else than the counting of linearly independent local operators classified 
by their ultraviolet scaling behaviour. 

3 Scaling Lee- Yang model with boundary 

The scaling Lee- Yang model with boundary is a combined bulk and boundary perturbation 
of the boundary version of the M.2.5 Virasoro minimal model which was investigated in 
detail in [16]. The conformal field theory has central charge c = —22/5 and the Virasoro 
algebra has two irreducible representations Vh with highest weight h = h^i = and 
h = hi t 2 = —1/5 [TZ]. We define the truncated characters of these representations by 



aF + f3f = {aF n + {3F n } 




(2.8) 



X 



00 




n=0 
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where d{n) gives the degeneracy of the level n descendents. They can be represented as 
the following fermionic sums [18] 



Xl,l(<?) ^ ^ _ q2+5n^l _ g3+5n) ' Xl,2(?) (1 _ gl+5n) ^ _ g4+5n) ^-1) 

Boundary conformal field theory was developed in \19\ EHJ EE] and the interested reader 
is referred to them for details. Applying the formalism to the conformal Lee- Yang model it 
can be seen that there are two conformally invariant boundary conditions. On one of them, 
denoted by 1 in [16], the spectrum of boundary fields is given by the vacuum representation 
Vo of the Virasoro algebra, and therefore it does not have any relevant boundary fields - 
thus can have no boundary perturbation either. In the other case, denoted $ in [T6] . 
the spectrum of boundary fields is given by the direct sum Vo © V-1/5 and therefore it 
has a nontrivial relevant boundary field (p with scaling dimension —1/5 and the general 
perturbed boundary conformal field theory action can be written aa^l 



Ax^h) = A^ + fj. I dy I dxc/)(x, y) + \i B \ dy<p(y) 



where Aj> denotes the action for .M(2/5) with the $ boundary condition imposed at x = 0, 
and \x and \ib denote the bulk and boundary couplings respectively. There is a unique 
nontrivial relevant bulk perturbation given by the spinless field with scaling dimensions 
h — h = — 1/5. The action of is similar, but the last term on the right hand side is 
missing. For fi > the bulk behaviour is described by an integrable massive theory having 
only a single particle with mass m with the following S matrix [22]: 



s i n h ( I + HE£) 

(x)= )l 2 { [x] = {x){l-x) 
V ' sinh (| - if) ; 



The pole at 6 = corresponds to the "<£> 3 property", i.e. the particle appears as a bound 
state of itself. The minimal bulk two-particle form factor only has a zero at 6 = and a 
pole at 6 = in the strip < ^sm(9) < n and is of the form [23] : 

f{6) = ^^v{m - B)v{-m + 9) , y = e e + e~ e (3.2) 



where 



. . dt • et sinh % sinh | sinh 
v{9) = exp <^ 2 / — e l ~ 



t sinh 2 1 



3 The field L_iip = d y (p, albeit relevant, does not induce any nontrivial perturbation because it is a 
total derivative. 
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and it satisfies 



sinh 9 



sinh 9 — i sin | 



f{0)f{0 + i*) = 

tt\ cosh# + l/2 
cosh 9+1 



f{o + il)f(e-i 



1(0) 



and consequently also f{9) ~ 0(1) for large 9. 

In the boundary theory with the perturbed <3> boundary, the reflection amplitude of the 
particle depends on the strength of the coupling constant of the boundary field as [16J 



R(9)* = Ro(9)R(b,0) 
with the one-particle coupling 



6+1 



where the relation between the dimensionless bootstrap parameter b and the dimensionful 
Lagrangian parameters /x and /is is known explicitely |24| . while in the case of the 1 
boundary the reflection amplitude is the parameter independent expression 



R(9)i 



and the one-particle coupling is 



9l = -z2(3) 1 / 4 (2 - V3) 



1/2 



3.1 Counting towers for the boundary condition 1 
3.1.1 The recursive equations for the form factors 

The 1PFF function corresponding to the parameter free reflection factor R{9)\ is 

n(9) = ismh9 u(9), 

where 



u{9) = exp 



00 dt 
T 



sinh | 



t 






-COS 


[(t-» 


n 









sinh — + sinh ~ — sinh 
sinh 2 £ 



t -i 



and asymptotically behaves as r t ~ e 2e when 9 — > oo. Following the general ideas in 2.2.1 
we introduce the Ansatz 

F n (0 1; . . . A) = 4"iWyi, • • • , y») II — II m ~ Bi) !^ i + ei) (3.3) 



!/* + % 



where 



H„ 



2^(0) 

Then one finds the following recursion equations for the Q n -s [5]: 
V: Q 2 (y + ,y-) = Q 1 (y) 

n 

Q n+2 {y+,y-,yi, ... ,y n ) = Ufa + 2/0 Qn+i(y,yi, ■ ■ -,y n ) > n > °; 



(3.4) 



8=1 



/C: Q 2 (-y,y) = 

Qn+2{-y,y,yi,---,yn) = Pn(y\yi,---y n )Qn(yi,---,yn) , n>0; 

B : Qxfa) = 

Q n+1 (0,y 1 ,...y n )=B n (y 1 ,...,y n )Q n (y 1 ,...,y n ) , n > 0; 



where 



2{y+-y-) 



Y[(yi - y-)(yi + y+) - Y[(y% + y-)(v% - y^ 



i=i 



1 / n n \ 



and 



y+ 
y- 



UJZ + uj 1 z 1 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



u = e s 



with the auxiliary variable z defined as a solution of y = z + z~ x (i.e. writing y = 2 cosh^ 
we obtain z = e 9 ). The symmetry of the expressions (|3.5l 13.61) in y± ensures that the 
resulting relations do not depend on which of the two possible solutions of the relation 
y = z + z^ 1 is chosen (switching from z = e 9 to z = e~ 9 only interchanges y + with y_). 

Note also that Q n+ 2(y+, y~, yi, . . . , y n ) only depends on y rather than separately on 
y± because due to its symmetry in all variables it can only depend on the combinations 



y- = y and y+y_ 



3. Furthermore it is easy to see that P n is a polynomial 



in all of its variables (since the expression inside the parentheses vanishes if y + = y_ and 
is therefore divisible by y + — y_), and using the previous argument again it only depends 
on y. However, elimination of y± makes the formulae rather cumbersome and therefore we 
prefer to keep them in our equations. 

It is also interesting to note that the boundary kinematical recursion is actually "almost" 
entirely contained in the bulk recursion. Using the bulk kinematical recursion with y = 
results in 

Q n+2 (0, 0, s/i, ... , y n ) = P n (0\yi ■ ■ ■ y n )Q n (yi, ■■■,y n ) 
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while using the boundary kinematical recursion twice leads to 



Q n+2 (0, 0, y h . . . , y n ) = B n (0, y x ... y n )B n {y x . . . y n )Qn(yi, ■■■,Vn) 

and therefore it must be true that 

Pn(0\yi ■ ■ ■ Vn) = B n (0, y x ... y n )B n (y 1 ...y n ) (3.9) 

which is indeed satisfied by the expression (13.61 13.71) . This means that the "square" of 
the boundary kinematical recursion is contained in the bulk one, as already discussed in 
[5]. The only independent information that the boundary kinematical singularity axiom 
carries is the sign of the one-particle coupling g: all the other axioms depend only on the 
one-particle reflection factor R, which in turn contains only g 2 using eqn. (12. ID . 

Indeed it was shown by Dorey, Tateo and Watts in [25] that although the fundamental 
reflection factors of the boundary condition 1 and $(6 = 0) are identical, these boundary 
conditions are physically different. Their one-particle couplings differ by a sign 

g x = -i2 (3) 1/4 (2 - V3) 1/2 = -g Hb=0) 

and, as a result the whole bootstrap structure changes: the 1 boundary condition has no 
boundary excited state, while the $(6 = 0) does. This is also reflected in the different 
expression for their 1PFF function r as discussed after eqn. (13.101) . We note that while the 
sign of g is not manifest in the reflection factor itself, it affects many physical quantities 
besides the spectrum such as finite size corrections to the energy levels [16, 26J and one- 
point functions of the bulk fields [3 127] . 

3.1.2 Asymptotic behaviour of form factors and simple towers 

It can be easily deduced that the asymptotic behaviour of the Ansatz ( 13.31) is given by 

\F n {6 1 + e,... i e n + 6)\ ~ e x ^ where x n = deg Q n - " ( " ~ 3) 

9 —> +oo Z 

where deg P denotes the total degree of the multivariable polynomial P defined by 

P(Xy u ...,Xy n ) ~ A dcgP as A -> oo 

Note that 

deg£>(Q n+ i) = degB(Q n+1 ) = degQ„ +1 - n + 1 , deg/C(Q n+2 ) = degQ n+2 - 2n + 1 

As a result, the solutions of the recursion relation can be expanded in the basis of simple 
towers which are defined to be towers with the property 

Q n = : n<n 
deg Q n+1 = deg Q n + n - 1 : n>n 
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for some n G N. 

For a simple tower the exponent x n describing the asymptotic growth of the form factor 
tower is independent of the level n. Therefore we can use the naive dimension to introduce 
a grading on the space of towers: the subspace of towers of grade x is spanned by the 
simple towers with naive scaling dimension x. For a simple tower given by the polynomials 
{Q n } neN the scaling dimension can be expressed as 

x = deg Q n 

3.1.3 The recursion kernel and the number of independent towers 

Here we carry over an approach developed by Koubek [TOl CLE] to count solutions of the bulk 
form factor equations to the boundary case. The idea is to classify the so-called 'kernels' 
of the recursion relations (12.61) first and use these to get the dimensions of the spaces of 
simple towers. 

The kernel can be defined as polynomials (at a given level n) that are taken to zero 
by the appropriate recursion relation in ( 12.61) . The kernel of each recursion at level n can 
be generated by multiplying an elementary kernel polynomial with an arbitrary symmetric 
polynomial of n variables. The generating polynomials are 

£ : K n(Vli ■ ■ ■ , Vn) = IJ (Vi + Vj) 

l<i<j<n 

V : K%( yi ,...,y n )= J] (y-+yiyj + y--3) 



l<i<j<n 
n 



B : K»( yi ,...,y n ) = H 



Vi 



i=l 

The common kernel of the three recursions at level n is generated by the least common 
multiple of these polynomials which is just equal to their product: 

K n (yi, ...,y n ) = K*(y 1} . . . , y n )K^( yi , y n )K°(y 1} ...,y n ) 

with degree 

n(3n-l) 
deg K n = 

and a basis of the kernel polynomials at level n is given by 

aff . . . a [ ^K n , < ki < k 2 . . . < h < n 

where are the elementary symmetric polynomials of n variables and degree k defined 
by the generating relation 



i=l k=0 
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The number of linearly independent towers can be counted using the observation that 
every independent tower starts from a kernel polynomial, because each new tower appears 
at some level (number of particles) n as an ambiguity of the solution of the recursion 
relations. This is true even for the towers starting at the lowest level n = 1 if we note that 
the kernel at level 1 is generated by 

#1(2/1) = 2/1 

The naive scaling dimension of the simple tower starting from the polynomial . . . a^K n 
is 

4:U = deg . . . o%>K n - = n(n + 1) + + + *, 

Now we can write down the generating function X(q) defined in (12.81) easily: 



00 00 

n(n+l)+m 



x (q) = E E p( m \ n )i r ' 



n=l m=0 



where p(m\n) denotes the number of partitions of the number m such that none of the 
summands is greater than n. Using the identity 



00 n 1 

x>Mn)<r=n— 



m=0 i=l 

we obtain 

00 „n(n+l) 

*<«) = E-I — 

n=1 Y[ 1 - q* 

i=i 

and with the help of a Rogers-Ramanujan identity (following [XT] ) 



x 1 

X ^ = E ( 1 _ q 2 +5 n )(1 _ q 3 +5 n ) ~ 1 = XlM 
n=l ^ ' 



using (I3.ip . Therefore the result of the counting of the form factor towers, graded by their 
naive scaling dimension, exactly matches the conformal spectrum of the operators localized 
on the 1 boundary, which is given by the vacuum representation of the Virasoro algebra 
(with highest weight h\,\ = 0). The subtraction —1 corresponds to the identity operator 
which only has trivial (vanishing) form factors. 

Note that in claiming the agreement we tacitly supposed that the naive scaling dimen- 
sion x of the towers can be identified directly with the conformal weights of the scaling 
operators, i.e. that there are no anomalous dimensions. This was checked for the form 
factor tower with the lowest possible dimension x = 2 in [5], where performing the spectral 
sum we obtained that the ultraviolet dimension of the corresponding operator was indeed 
A = 2. 
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3.1.4 Some explicit solutions 



Finally we give the explicit forms of the lowest levels of some form factor towers. The 
n = 1 generating kernel polynomial is K\ — a\ , and all kernel polynomials can simply be 
written as cr™. The first three towers starting at level 1 are 

x = 2 : Ql = (7i , Ql = a l , Qj = o\ , Q\ = a\{o 2 + 3) ... 



x = 3 
x = 4 



r>9T _ „ n T 

Vn — ^lVn 
-)d 2 T 2nT 



The upper index T shows that the tower corresponds to the (off-critical version of the) 
boundary stress-energy operator T = L_ 2 I (already found in [5]), while the other towers 
are its first and second derivatives respectively. In fact, there is a single tower for each value 
of x less than 6 which corresponds to the fact that up to this level (due to the presence 
of null- vectors) the conformal vacuum representation of the c = —22/5 Virasoro algebra 
contains only the vectors L™ 1 L_2l {n = 0, 1, 2, 3). At x = 6 the conformal representation 
contains another linearly independent vector, and indeed a new tower appears, which starts 
from the level 2 generating kernel polynomial: 

Qf = a i a 2 {<?l -o- 2 -3) 

Qf^ = 0i (a\a 2 - cr 3 ) (a\ - cr 2 — 3) 

Q { t = °"i (°2 + 3) ((o"i<t 2 - 03) (<jI - cr 2 - 3) + 0-104) 

The counting argument in the last subsection ensures that the new towers always appear 
in such a way that the space of simple towers graded by x matches the space of linearly 
independent vectors in the vacuum representation. 



3.2 Counting towers for the boundary condition $ 

In this case the 1PFF is [5] 

isinh# ,„ s , 

r M = ^rrr M W ( 3 - 10 ) 

(sinh 9 — i sin w( 6 J )(sinh 9 — i sin n „ ' ) 

The main difference from the I case is the presence of the denominators, which correspond 
to boundary excited states. Note that at b = the reflection factor R{9)<& is identical to 
R(9)i, but the corresponding lPFFs are different. This is related to the different interpre- 
tation of the pole structure of the R matrices discussed in detail in [25]: for the boundary 
condition 1 the boundary excited state is absent and the corresponding pole at 9 = i| is 
explained by a Coleman-Thun diagram. In this case we take the following Ansatz 

x? (a a\ u n ( 77 r $(^) TT [(& ~ + ft?) /qua 
F n (9 1 , ...,6 n ) = H n Q n (y u ...,y n ) [[ [ [ (3.11) 

i Vi i<j Vi + Vj 



4 We omit the upper index n of specifying the number of variables since it is always the same as 
the number of particles, i.e. the level of the form factor tower. 
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The recursion relations in this case are [5] 
V: Q 2 (y + ,y_) = (y 2 -^ 3 )Q 1 (y) 

n 

Q n +2(y+,y-,yi,---,yn) = (y 2 -f3- 3 )Y{(y + yi)Qn+i(y,yu---,yn) , n > 0; 

i=l 

/C: Q 2 (-j/,j/) = 

Q n +2(-y,y,yi, ...,y n ) = (y 2 - l3-i){.y 2 - (3f)P n (y\yi, ■ ■■y n )Q n (yi, ■■■,y n ) , n>0; 
B : Qi(0) = 

Qn+lCOjI/l* • • -3/n) = P-lPlBnfa, . . -yn)Qn{yU ■ ■ ■ ,Vn) , n > 0; 

where P n and S„ are given by eqns. (13.61 13.71 13.8ft and 

(3 k (b) = 2 cos ^(b + k) 
o 

The bulk and boundary kinematical recursions again satisfy the appropriate generalization 
of the compatibility relation (13.91) . 

For the degrees of the recursions we have 

deg£>(Q n+ i) = deg B(Q n+1 ) = degQ n+1 - n- 1 , deg/C(Q n+2 ) = degQn+2 - 2n- 3 

and so simple towers can now be defined by the property 

Q n = : n <n 
deg Q n+ i = deg Q n + n + 1 : n > n 

Due to the additional factors in the denominator the asymptotic behaviour of r${6) for 
9 — > oo is 0(1), and so the naive scaling dimension of a simple tower corresponding to the 
polynomials {Q„} neN is 

x = deg Q n - n ( n ^ L 

The kernel is generated by the same polynomials as before, but now the naive scaling 
dimension corresponding to . . . cr? K n is 

4:U = deg^ • --^Kn - = n{n-l) + k 1 + ... + k l 

and result we obtain 

x x 00 qn(n-l) 00 qn{n+l) 00 ^n 2 



m = EE^H^ n(n - 1)+m = E^ = + £ 



if 



n=l m=0 

i=l i=l i=l 



1 JJ 1 - q i n=1 Y[ 1 - q i n=1 Y[ 1 - g< 

8=1 

+5n\ n — q 1+5n )(l 



(1 - Q 2 + 5 ")(l - 3 + 5 ™) ^ (1 - Q 1 + 5 ")(l - Q 4 + 5 ™) 
n=l v ' v ' n=l v ' v ' 
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using again a Rogers-Ramanujan identity following [TT]. This is consistent with the fact 
that in the conformal limit the boundary operator content of the <3> boundary condition is 
given by the direct sum of the vacuum module with highest weight hi t \ = and the only 
other nontrivial representation with highest weight hip = —1/5. 

Note that in this case to achieve the agreement we must suppose that the naive scaling 
dimensions of the towers corresponding to operators of the identity representation do not 
get any anomalous corrections, while those corresponding to operators in the model with 
highest weight hip = —1/5 get exactly the right contribution for their conformal weight 
to match with that predicted by conformal field theory. The true scaling weight of the 
operator can only be computed by examining appropriate two-point functions. This makes 
the operator identification very difficult for higher levels as shown below using the example 
of the boundary stress-energy tensor. 

For x = the operator identification is easy, since there is a unique tower with the 
required asymptotic properties: 

Q\ = a u Q% = ai(a 2 +3-i3 2 _ 3 ), 

Qt = d [01(02 + P-i)(v2 + Pi)- (^2 + 3)(a l( x 2 - o 3 )] (3.12) 

In [5] the two-point function of this tower was compared to the perturbed conformal field 
theory prediction and we found that its identification with the operator ip is justified (this 
also settles the case for all the derivatives of <p). We remark that the anomalous dimension 
of the correlator in this case turns out to be 

2A = — ^2/i=-- 
5 T 5 

which is due to the fact that p itself appears in the pp operator product, and so the above 
result is in accordance with the discussion under 2.2.2 since h m[n = — 1/5 (cf. also [5] and 
also the analogous bulk situation discussed by Zamolodchikov in [23]). For x — 1 there is 
again a single operator 

Q d n = ^Ql 

whose scaling dimension can be known exactly in terms of the scaling dimensions of the 
lowest tower, since it is a derivative operator. However, we find two towers at x = 2. One 
is the second derivative 

Qn * = °\Ql (3-13) 

while the other one starts from the n = 2 generating kernel polynomial: 

Ql = ClOv, (of - o 2 - 3) , Ql = Ol (0102 - 3 ) (of - 2 - 3) + O1O3) 

which is consistent with the fact that both the Virasoro module h — 1/5 and the module h = 
have a single linearly independent vector at level two: L 2 _ x ip and T = L_ 2 II, respectively. 
The tower (13.131) can be directly identified with L^ip since the action of L_i is exactly 
given by the derivative, but the identification of T is more involved. In general we expect 
that 

T <xf + oup + pdip + i&y (3.14) 
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(where the proportionality sign means that the normalization of the operator must also 
be fixed). Note that T cannot mix with higher towers because that would spoil the short- 
distance behaviour. The identification of T is an open question, to which we return in the 
conclusions. 



4 Sinh-Gordon model with Dirichlet boundary condi- 
tions 

The sinh-Gordon theory in the bulk is defined by the Lagrangian 



m 



(cosh 6$ - 1) 



It can be considered as the analytic continuation of the sine-Gordon model for imaginary 
coupling (3 = ib. The S-matrix of the model is 



S(9) 



B 

1+ 2 



B 

T 



B 



B 



lb 2 



57T 



b 2 



In this case there is no self-fusion pole, and so bulk dynamical singularities are absent. The 
minimal bulk two-particle form factor belonging to this S-matrix is |28j 



/(0) = Afexp 



dx 



sin 



x 



x{m-9)\ sinh^f sinh(l - f)f sinhf 



2tt 



sinh x 



satisfying 



f(9)f(9 + i<ir) 



sinh 9 



sinh 9 + i sin ^ 



(4.1) 



and therefore f{9) ~ 0(1) as 9 — ► oo. 

Sinh-Gordon theory can be restricted to the negative half-line, but integrability is only 
maintained by imposing either the Dirichlet 

$(0,i) = $ D 

or the two parameter family of perturbed Neumann 



V B ($(0,t)) = Mo cosh ( -($(0,t) - $ ; 



boundary conditions [3J. The latter interpolates between the Neumann and the Dirichlet 
boundary conditions, since for M Q = we recover the Neumann, while for M — > oo the 
Dirichlet boundary condition with §q = $ . The reflection factor which depends on two 



continuous parameters can be written as [29 
R{9) = Rq{9)R{E,F,9) = 



1 B 
2 + 4 



1 - 



B 



E 
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and can be obtained as the analytic continuation of the first breather reflection factor in 
boundary sine-Gordon model which was calculated by Ghoshal in [30]. The relation of the 
bootstrap parameters E and F to the parameters of the Lagrangian is known both from 
a semiclassical calculation [2H [31] and also in an exact form in the perturbed boundary 
conformal field theory framework [32]. Imposing Dirichlet boundary condition instead of 
the general one corresponds to removing the F dependent factor from R{9). 



4.1 Recursion relations for the form factors with Dirichlet bound- 
ary conditions 

Boundary form factors for sinh-Gordon theory with Dirichlet boundary condition have 
already been investigated to some extent in [5], and many solutions were constructed by 
Castro-Alvaredo [33] at the self-dual point B — 1. Here we aim to classify all the solutions 
for general coupling. 

The 1PFF corresponding to Dirichlet boundary condition is [5] 



r E (9) 



sinh 9 



sinh 9 — i sin 7 



u(9,B) 



where 

u(9, B) = exp 
and 



, 2 1 dx 

x 



,171 



cos 



£-1 



cosh £ 



— fsinh h sinh(l ) — h sinh — 



B'X 



7r J sinh x \ 



2' 2 



71 



7= s (*-i: 



The asymptotics of the 1PFF is te ~ e e as 9 — > 00. The one-particle coupling to the 
boundary is given by 



9e 



2(l + cos^ + sin^f) cos 7 



sin 



tB 



1 — sin 7 

Writing the n-particle form factors in the general form 

r E {e i ) T ^f{9 l -9 3 )f{9 l + 9 ] ) 



F n {9 u ...,0 n ) = H n Q n ( yi , ...,y n )H H 



i<j 



Vi + Vj 



and choosing the ratio of the H n -s appropriately the recursion equations take the form: 

K: Q 2 (-y,y)=0 

Qn+2(-V,V,Vi,---,Vn) = (Z/ 2 -4cos 2 7)P n (y|y!, ...,y n )Q n (y 1 ,...,y n ) , n > 0; 
B : Qi(0) = 

Qn+i(0,yi,---,yn)=2cosyB n (yi,...,y n )Q n (y 1 ,...,y n ) , n > 0; (4.2) 
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where P n is given by eqns. (13.61 13.81) with to = e m 2 , 

B n ( yi ,...,y n ) = 4s /^ (f[ (^- 2sin ^) -]l(^ + 2sin ^) 

and the bulk and boundary kinematical recursions again satisfy the appropriate general- 
ization of the compatibility relation (13.91) . Note that for E = the right-hand side of the 
boundary kinematical relation is zero: 

B : Qn+i (0,?/i, ...,y n ) = 

and so one could dispense with the boundary kinematical recursion choosing an Ansatz of 
the form 



F n (9t,..., 9 n ) = H n Q n ( yi , ...,y n ) JJ^=o(^) fl 



f{e l -e J )f{e l + e ) 



Vi + Vj 

as in [5]. In this paper, however, we choose to treat the two cases together. 

4.2 Counting towers and the spectrum of local operators 

The recursion relations (I4.2f) have the following degrees: 

deg/C(Q n+2 ) = degQ n+2 - 2n - 1 , deg£(Q„ + i) = degQ n+1 - n 
and so simple towers are defined as 

Q n = : n <n 
deg Q n+1 = degQ n + n : n > n 

while the naive scaling dimension of a simple tower corresponding to the polynomials 
{Qn} n€N is given by 

x = deg Q n - - 

Due to the absence of bulk dynamical singularities, the generating kernel polynomial at 
level n is the product of the ones for the bulk and boundary kinematical recursions 

K'niVi, ■ ■ ■ ,Vn) = K^(y l7 . . .,y n )K^(y u ...,y n ) 

with degree 

deg< = ^ 

and the naive scaling dimension of the simple tower starting from the polynomial . . . a£ K' n 
is 

= deg 4? • • • - =n + k 1 + ... + kl 
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The generating function (12.80 takes the form 

OO OO OO OO 

x ^ = E E K™w<r m = E ^ n )<? n = n 

n=l m=0 n=l j'=l 

where jp(^) is simply the number of all partitions of n. 

On the other hand, the independent boundary operators for the Dirichlet boundary 
condition $(£, x = 0) = $o are simply given by the differential monomials 

d k x ^(t,x)...d^(t,x)\ x=Q 

which have scaling weight ki + . . . + ki. Therefore the number of operators with a given 
scaling weight n e N is indeed equal to the number of integer partitions of n. 



4.3 General boundary conditions 

When both E and F are nonzero, the 1PFF reads 

tef{9)= - . . ^ - . . - u(B,B) , 7=f(£-l) l' = ^{F-l) 
(sinn — « sm 7j (smh t7 — zsm7'J 2 2 

with asymptotics r Si 7 ~ 0(1) when 9 — > oo, and therefore the scaling dimension of a simple 
tower is 

x = deg Q n - n ^ 2 + - 

and so each tower starting from an elementary kernel polynomial K' n has naive scaling 
dimension x = 0. This results in infinitely many towers corresponding to any integer value 
of the naive scaling dimension! 

This is in fact not so surprising if we consider the expected spectrum of scaling operators. 
These are of the form 

<9 t fcl $(t,x) ...d^(t,x)e Q ^ x) 

x=0 

(the ^-derivatives of the field can be expressed with exponential operators using the bound- 
ary condition). These can be thought to be organized into families of descendents of the 
exponential operators 

e a>S>(t,x=0) 

the descendent level given by ki + . . . + ki, so the number of descendents at some level n 
is the number of partitions of the integer n. However, the exponential operator has naive 
(classical) scaling dimension zero: its scaling weight is a fully quantum effect. Although 
at first it seems that there exists a continuum of such operators, recall that they can be 
expressed in terms of powers of the field: 



k=l 
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which shows that they depend only on countably many independent operators, the powers 
$ fc (which are not scaling fields in themselves, their correlation function involves loga- 
rithms) so one expects only a countable infinity of form factor towers with naive scaling 
dimension 0, which is indeed what we found. 

In this paper we cannot go further with the identification between the form factor towers 
and the local operators because some further tools are necessary in order to classify the 
towers of zero naive scaling dimension arising from the kernel solutions; see the conclusions 
in the next section for details. We only wish to note that our results in the free field limit 
are identical to those in [5j, and that the whole situation above is very similar to the case 
of the bulk sinh-Gordon model [9|. 

5 Conclusions and open problems 

In this paper we performed the classification and counting of the solutions of the bound- 
ary form factor bootstrap based on their ultraviolet behaviour. We showed that in the 
boundary version of the scaling Lee- Yang model and in the boundary sinh-Gordon model 
the results are in perfect agreement with the expectations based on boundary conformal 
field theory in the first case and on the Lagrangian formulation in the second. This gives 
an additional and crucial piece of evidence for the consistency of the bootstrap axioms 
developed in [5]. 

Our discussion, however, has made several open problems manifest, that are worthwhile 
to be pursued in the future. 

5.1 The issue of anomalous dimensions 

The first is the question of anomalous dimensions. It was shown (using the example of the 
boundary operator ip and its descendents in the Lee- Yang case) that the naive dimension 
obtained from the asymptotic growth of the form factors of some tower is not necessarily 
identical to the exact ultraviolet weight of the corresponding operator. In fact the counting 
of towers depends very much on the additional assumption that anomalous contributions 
to the scaling dimensions do not mess up the spectrum_|, and so the space of simple towers 
with a given value of the naive scaling dimension x can be brought into correspondence 
with the space of local operators of a given value for the exact ultraviolet weight h. 

This issue is also well-known in the bulk and there the computation the ultraviolet 
weight can be tackled by two different approaches. One of them uses the cumulant ex- 
pansion of the logarithm of the two-point function [34] (for a very nice discussion see also 
[35]). The spectral expansion (12.41) can be written in the following form for the Euclidean 

5 It could happen, at least in principle, that the grading of the operators is changed by the anomalous 
dimensions, so that the classification according to conformal descendents does not match the grading 
provided by the naive scaling dimension. 
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two-point function 



00 1 poo poo poo 

p(mr) = / d0 i d9 2 ... d9 n e- mT ^ coshe *f n (9 1 ,...,9 r , 

n _Q J —00 J —00 J —OO 



1 

(47) 



f = F F+ 



where we used that the functions f n are symmetric and even in all the rapidity variables 
due to the symmetry properties of the form factors (note that in a unitary theory they are 
also positive semidefinite, but that is not necessary for the following argument). Supposing 
that the leading n = term in the short-distance limit is a constant and normalizing it to 
1 we can write a similar expansion for the logarithm of the correlation function 



00 -1 POO POO POO 

logp(mr) = J2 - / Mi d9 2 ... d9 n e- mT ^ cosh0 <h n (9 U 
„_i J —00 J —00 J —00 

where the h n are the cumulants of the functions f n defined recursively by 



On) 



Wx) = h{9 1 ) , f 2 (ei,6 2 ) = h 2 (6 1 ,6 2 ) + h 1 (9 1 )h 1 (9 2 ) 

f 3 (0i, 2 , 9 3 ) = h 3 (6 u 6 2 , 9 3 ) + h 1 (6 1 )h 2 (6 2 , 6 3 ) + h^h^t, 9 3 ) + h^h^, 6 2 ) 
+h 1 (9 1 )h 1 (9 2 )h l (9 3 ) 



Defining 

h n (9 1 ,...,9 n )= \imh n (9 1 + X,...,9 n + \) (5.1) 

A^oo 

it is easy to see that the functions h n depend only on the differences of the rapidities 
(for spinless operators in the bulk this is already true for f n and therefore also h n and in 
that case h n = h n ). Following the same argument as in the bulk case we arrive at the 
representation 

00 1 f°° f°° 
A = J2~ d9 2 ... d9 n h n (0,9 2 ,...,9 n ) (5.2) 

n=l n ' 

for the short-distance exponent defined in d 2 . 5 ft . 

The main problem with this approach is that the limit (15.11) exists only for boundary 
operators for which the naive scaling dimension is zero, and therefore it is of very limited 
use (in the bulk it obviously exists for any spinless operator, but a similar problem appears 
when trying to apply this approach to operators with non- vanishing Lorentz spin). An 
example to which the cumulant expansion (15 .2p can be applied is given by the operator tp 
in the scaling Lee- Yang model with $ boundary condition, but in that case we already have 
a much more detailed comparison with the conformal prediction directly via the two-point 
function, which was carried out in [5]. Therefore it is an interesting problem to develop 
some method to extract the ultraviolet dimension directly using only the spectral densities 
f n , generalizing the above argument. 
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The other approach in the bulk is to use a sum rule such as the so-called A-theorem 
developed in [36]. Unfortunately, the arguments used in the bulk case are not directly ap- 
plicable to boundary theories, because the boundary component of the energy-momentum 
tensor does not obey any conservation law in itself, but it remains to be seen whether there 
exists some other way of deriving a sum rule. 

A further possibility along this line of thought would be to examine the operator prod- 
uct of T (the operator that corresponds to the conformal boundary Virasoro current off 
criticality) with the boundary operator O in question using the spectral representation for 
their two-point function, since it is known from boundary conformal field theory that they 
have an operator product expansion of the form 



but to perform this we would need to identify the form factors corresponding to T; the 
discussion below just highlights the difficulties related to operator identification. 

5.2 Operator identification 

It is also an open problem how to identify individual operators with their corresponding 
towers. In some cases this is straightforward, like for the operator with the lowest conformal 
weight (an excellent example is again the operator ip in the scaling Lee- Yang model with <E> 
boundary condition, or the boundary stress energy tensor T in the case of the 1 boundary 
condition j5j). It is also rather easy to identify form factors of derivative operators since 



(once the form factors of O are known). For other operators, however, the identification 
is not at all straightforward. Take the example of the boundary stress energy tensor T 
for the $ boundary condition in the scaling Lee- Yang model (subsection 3.2) where the 
new tower at x = 2, corresponding to the appearance of T, mixes non-trivially with other 
towers according to (13.140 . It is not easy to disentangle the mixing and so produce a 
direct identification. There are certainly some approaches to try which we leave to future 
investigations. 

One possibility is to calculate the mutual two-point correlation functions of the towers 
T, <f, dp, d 2 ip and compare it numerically to the conformal perturbation theory prediction 
for the operators T, <p, L_\(p, L 2 _ x p> (in fact only the correlators involving T are necessary). 
Using that, one can then fit the mixing (and normalization) coefficients numerically. 

Another approach to operator identification in general is to use a boundary extension of 
the approach developed by Delfino and Niccoli for the bulk scaling Lee- Yang model in [37] 
who used the fact that part of the conformal descendents can be generated using derivation 
{L-x) and the charges which remain conserved off-criticality. This is very efficient for low 
levels: it makes possible the precise identification of towers with the appropriate operators 



T(t)O(0) 



h o O(0) 

T 2 



+ less singular terms 
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up to descendent level 7. It remains to be seen whether this method can be extended to 
the boundary case. 

In general, however, it is also difficult to identify towers corresponding to the primary 
fields. As we have shown in subsection 4.3, this is a very complicated issue in the boundary 
sinh-Gordon model with general boundary conditions, where an infinity of x = towers 
must be matched against the spectrum of exponential fields. In the bulk case form factors 
of primary fields satisfy a factorization property 0, [34], [38], for which a general argument 
was given in [36]. Factorization makes possible the identification of towers corresponding 
to primary fields, and can also be extended to descendents. It remains to be seen whether 
such clustering property can be extended to the boundary form factor bootstrap: at first 
sight, it seems more promising to take an approach which relates clustering to the operator 
product expansion as in the paper [39] by Balog and Weisz instead of one appealing to 
conformal holomorphic factorization (as done by Cardy et al. in [36]), which is broken by 
the presence of the boundary condition. In this regard it is also interesting to examine the 
relation between the bulk and boundary form factor bootstrap, to which we now turn. 

5.3 Connection to the bulk form factor bootstrap 

Isolating the leading asymptotic coefficient of the n-particle form factor F® in the large 
rapidity limit 

F n °(9 1 + A, . . . , n + A) ~ F® (0 U . . . , 9 n )e xX + subleading terms 

it is obvious that F®(0i, . . . , 9 n ) is a function of rapidity differences only, and it is also easy 
to verify that it satisfies the bulk form factor axioms as a consequence of F® satisfying the 
boundary form factor axioms. Therefore F®(9i, . . . , 6 n ) can identified with the n-particle 
form factor f®(9i, . . . , 9 n ) of some bulk operator O. As an example, using the asymptotic 
behaviours of the 2PFF function / in (13 .21) and of the 1PFF function r$ in (13 . 1 Oh . the form 
factor tower F v (13.120 in the Lee- Yang model with <3> boundary condition in this limit 
gives rise (up to normalization) to the bulk form factor tower of the trace of stress energy 
tensor 9 found by Zamolodchikov [23]. As an application of this correspondence it is easy 
to see that the bulk cumulant expansion for the anomalous dimension of 9 in [34] is just 
twice the expression (15.21) evaluated for the boundary field (p. This is consistent with the 
fact that in the bulk A = — 2/5 while for the boundary case = —1/5. 

This correspondence may be related in some way to the bulk-boundary operator product 
expansion in boundary conformal field theory [20], and it is possible that (extending the 
example above) the exact scaling dimension of operators can be identified this way using 
known results for the bulk form factors. It could also explain why the detailed structure of 
the counting of form factor towers (kernels, naive scaling dimensions, character identities 
to use etc.) is so similar to the bulk case as it was treated in [9] HOj CLE]. 
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